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We consider the conditions which must be satisfied for a Majorana RH sneu- 
trino, a massive right-handed (RH) sneutrino associated with the see-saw mech- 
Q_i' anism of Majorana neutrino masses, to play the role of the curvaton. Planck- 

• scale suppressed non-renormalizable terms in the RH neutrino superpotential 

must be eliminated to a high order if the RH sneutrino curvaton is to dom- 

X 

5_i ■ inate the energy density before it decays, which can be achieved via an R- 



symmetry which is broken to R-parity by the Giudice-Maseiro mechanism. In 
order to evade thermalization of the curvaton condensate, one RH neutrino 
mass eigenstate must have small Yukawa couplings, corresponding to a lightest 
neutrino mass m Ul < 10~ 3 eV. A time-dependent lepton asymmetry will be 
induced in the RH sneutrino condensate by the Affleck-Dine mechanism driven 
by SUSY breaking B-terms associated with the RH neutrino masses. Requiring 
that the resulting baryon asymmetry and isocurvature perturbations are accept- 
ably small imposes an upper bound on the RH neutrino mass. We show that a 
scenario consistent with all constraints is obtained for a RH neutrino mass in the 
range 10 2 — 10 GeV when the RH sneutrino decay temperature is greater than 
100 GeV and lightest neutrino mass m Ul ~ 10~ 3 eV. Larger RH neutrino masses 
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are possible for smaller m Ul . The resulting scenario is generally consistent with 
a solution of the cosmic string problem of D-term inflation. 
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1 Introduction 



The possibility that the density perturbations responsible for structure formation come 
from quantum fluctuations of a scalar field other than the inflaton has recently at- 
tracted considerable interest [U El El ■ The scalar field, known as the curvaton [3], 
is effectively massless during inflation, but begins coherent oscillations after inflation 
and comes to dominate the energy density of the Universe before it decays. The initial 
isocurvature perturbations of the curvaton are then converted to adiabatic energy den- 
sity perturbations. A number of curvaton candidates have been proposed and analyses 
of curvaton dynamics performed |U El El El El El UHj- 

An important question is whether there is a natural curvaton candidate in the 
context of the minimal supersymmetric (SUSY) Standard Model (MSSM) and its ex- 
tensions. SUSY is advantageous as it allows us to understand the flatness of the 
curvaton potential with respect to quantum corrections. A specific motivaton for a 
SUSY curvaton is provided by D-term hybrid inflation jllj . In D-term inflation with 
naturally large values of the superpotential coupling, the energy density during infla- 
tion must be lowered in order to evade the problem of cosmic microwave background 
(CMB) perturbations due to cosmic strings This in turn requires an alternative 
source of the cosmological energy density perturbations, such as a curvaton. 

Some general conditions for a SUSY flat direction scalar to play the role of a 
curvaton were considered in jB] . Similar conditions for the specific case of right-handed 
(RH) sneutrinos (the scalar superpartners of the RH neutrinos) were discussed in 
It was shown that small couplings of the curvaton to the MSSM fields are required in 
order to evade thermalization of the curvaton condensate by scattering from particles 
in the thermal background and in order to have a low enough decay temperature for 
the curvaton condensate to dominate the energy density before it decays. Curvaton 
domination also requires the elimination of Planck-scale suppressed non-renormalizable 
terms (NRT) to a high order in the flat direction superpotential. 

The case of a coherently oscillating MSSM flat direction condensate formed from 
a linear combinations of squark, slepton and Higgs fields was considered in j^j. It was 
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shown that as a result of their interaction with the fields in the background radiation, 
the condensates are either thermalized or fail to dominate the energy density before 
they decay. MSSM flat direction scalars can therefore play the role of a curvaton only 
if one assumes that the inflaton can primarily decay to 'hidden radiation' consisting 
of hypothetical light hidden sector fields [H]. In the following we will consider the 
conventional case of inflaton decay to observable radiation. 

An alternative curvaton candidate, in the context of the MSSM extended to acco- 
modate neutrino masses, is a RH sneutrino. The curvaton condensate is most likely to 
survive in the thermal enviroment of the early Universe and to dominate the energy 
density before it decays if it is relatively weakly coupled to MSSM fields. This makes a 
RH sneutrino a promising candidate for a SUSY curvaton. This possibility was noted 
in |H] and discussed in EI| , which also considered the possibility of leptogenesis via 
RH sneutrino decay and the associated isocurvature fluctuations. In [7] the conditions 
for a RH sneutrino to serve as a curvaton were considered in some detail. 

The results of [7j showed that a Dirac RH sneutrino (a RH sneutrino associated 
with Dirac neutrino masses, corresponding to RH neutrino masses Mn = 0), could 
serve as a curvaton if the amplitude of the RH sneutrino at the onset of its coherent 
oscillations satisfied N osc > 10~ 5 M, where M = M P i/\/&k with Mpi the Planck mass. 
This requires that Planck-scale suppressed superpotential NRT cx N n are suppressed 
for n < 6 [7]. The weak Yukawa coupling of the Dirac RH sneutrino ensures that 
there is no danger of condensate thermalization. However, the decay temperature 
of the RH sneutrino condensate was found to be around 10 MeV, which introduces 
two problems [JJ. The decay of the condensate RH sneutrinos will produce a large 
number of lightest SUSY particles (LSPs), which will overdose the Universe unless 
they decay sufficiently rapidly. In order that their decay products do not destroy 
the light elements produced in nucleosynthesis, the RH sneutrinos and LSPs must 
decay before nucleosynthesis occurs at T nuc « 1 MeV. This requires new R-parity 
violating interactions in order to have sufficiently rapid LSP decay. (This also has the 
disadvantage of eliminating the possibility of LSP dark matter.) The second problem 
is that the baryon asymmetry must be produced during the decay or after the decay 
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of the condensate RH sneutrinos, since any pre-existing baryon asymmetry will not 
be correlated with the entropy fluctuations produced by curvaton decay, resulting 
in unacceptably large baryon isocurvature fluctuations [5]. Therefore new baryon 
number violating interactions must be introduced in order that the RH sneutrino 
decay can directly produce the baryon asymmetry. Whilst no definitive analysis of 
these possibilities has yet been undertaken, the Dirac RH sneutrino curvaton appears 
to be disfavoured. 

Majorana RH sneutrinos (RH sneutrinos associated with the see-saw mechanism for 
Majorana neutrino masses [Hj, corresponding to the case ^ 0) can also dominate 
the energy density of the Universe before they decay. However, because of the larger 
Yukawa couplings in the see-saw mechanism for neutrino masses, the RH sneutrino 
condensate decays earlier, requiring a larger initial amplitude, N osc > 10~ 3 M, and 
correspondingly greater suppression of the Planck-suppressed superpotential NRT (all 
superpotential terms oc N n with n < 11 7\). In addition, the Majorana RH sneutrino 
curvaton was shown to have a potential problem with condensate thermalization. Due 
to its larger Yukawa couplings, the interaction of the Majorana RH sneutrino curvaton 
condensate with the MSSM fields in the thermal background is much stronger than in 
the case of the Dirac RH sneutrino. As a result, any contribution to the effective light 
neutrino mass matrix resulting from couplings to the curvaton's fermionic partner must 
be less than 10~ 3 eV in order to avoid thermalization of the condensate [7J. Thus if the 
masses of the neutrinos are characterized by the mass squared splittings suggested by 
solar and atmospheric neutrino oscillations, m v ~ 0.01 — 0.1 eV the RH sneutrino 
condensate would be thermalized before it could dominate the energy density. This 
will be the case if there is no correlation between the RH neutrino mass matrix and 
the neutrino Yukawa coupling matrix, in the sense that diagonalization of the RH 
neutrino mass matrix does not result in approximate diagonalization of the neutrino 
Yukawa coupling matrix. In this case the RH sneutrino curvaton mass eigenstate 
would have Yukawa couplings to the Higgs and lepton doublets characterized by the 
largest neutrino mass eigenvalue, m u ^0.1 eV, ruling out the Majorana RH sneutrino 
as a curvaton. 



3 



However, it is possible for mass squared splittings and mixing angles associated 
with solar and atmosperic neutrino phenomenology to be accounted for by just two 
RH neutrino mass eigenstates. For example, this is the case in the single right-handed 
neutrino dominance (SRHND) model of neutrino masses A third right-handed 
neutrino and its superpartner could then have arbitrarily small Yukawa couplings to 
the MSSM fields, allowing its condensate to evade thermalization. In this case the 
Majorana RH sneutrino becomes a viable curvaton candidate. The requirement that 
the lightest neutrino mass satisfies m vi ^ 10~ 3 eV may then be regarded as a prediction 
of the Majorana RH sneutrino curvaton model. 

In this paper we consider in detail the case of a Majorana RH sneutrino curvaton 
with a light neutrino mass eigenstate, which was not considered in [7j. An advantage 
of the Majorana RH sneutrino curvaton with respect to dangerous Planck- suppressed 
NRT is that the conventional R-symmetry of the MSSM (broken to R-parity by the 
fiH u Hd term) J7j, extended to include the RH neutrinos, will completely eliminate 
all RH neutrino superpotential terms except the mass term and Yukawa coupling 
terms responsible for neutrino masses, thus solving the problem of dangerous Planck- 
suppressed NRT. The R-symmetry must be broken to R-parity, but if this occurs 
via the Giudice-Masiero mechanism ^H], which generates only R-symmetry breaking 
effective superpotential terms with total R-charge equal to zero, then no dangerous RH 
neutrino superpotential NRT will arise. Operators eliminated by a global symmetry 
such as R-symmetry might be generated by quantum gravity wormhole effects |19j . 
However, this is dependent upon the short-distance structure of gravity [SU|. I n the 
following we will take the view that quantum gravity corrections are absent or highly 
suppressed. 

We will show that a time- dependent lepton asymmetry is induced in the RH sneu- 
trino condensate via a variant of the Affleck-Dine (AD) mechanism driven by the SUSY 
breaking B-terms associated with the RH neutrino masses. (This mechanism has also 
recently been discussed in the more general context of AD leptogenesis via a RH sneu- 
trino condensate in [SI].) Unlike the conventional AD mechanism associated with 
MSSM flat direction fields (SSj, L- violation from the B-terms remains unsuppressed 
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throughout the evolution of the condensate, resulting in a rapidly oscillating lepton 
asymmetry. If the RH sneutrino decay occurs before the electroweak phase transition, 
as is likely to be necessary in order to generate the baryon asymmetry, then any lepton 
asymmetry will be processed into a baryon asymmetry of the same order of magnitude 
by sphalerons, which must be less than the observed baryon asymmetry of the Uni- 
verse. The requirement that the lepton number to entropy ratio due to the decay of 
the RH sneutrino condensate satisfies ni/s'z 10~ 10 will be shown to impose a signifi- 
cant constraint on the Majorana RH sneutrino curvaton scenario. In particular, it will 
be shown that RH sneutrino leptogenesis via CP violating decays of the condensate 
RH sneutrinos ^3] can account for the baryon asymmetry only if the lightest neutrino 
mass is less than about 10~ 8 eV. This favours electroweak baryogenesis [231 as the 
origin of the baryon asymmetry in this scenario. 

The paper is organized as follows. In Section 2 we discuss the conditions for 
a Majorana RH sneutrino to dominate the energy density and to account for the 
observed energy density perturbations. We also discuss the possibility that a Majorana 
RH sneutrino could account for the energy density perturbations in D-term inflation. 
In Section 3 we discuss the constraints on the RH neutrino Yukawa couplings and 
neutrino masses coming from evasion of condensate thermalization. We also discuss 
the possibility that thermal effective masses could play a role in the evolution of 
the condensate. In Section 4 we consider the AD lepton asymmetry induced in the 
Majorana RH sneutrino condensate by the B-term AD mechanism and the associated 
constraints on the model. In Section 5 we discuss the range of model parameters 
allowed by the constraints. In Section 6 we present our conclusions. 

2 Curvaton Dominance and Density Perturbation 
Constraints 

We consider the simplest inflation scenario, in which the expansion rate during in- 
flation, Hj, is constant and inflation is followed by a period of coherent inflaton os- 
cillations. (We refer to this period as inflaton matter domination, IMD). Inflaton 
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oscillations decay away completely at the reheating temperature, Tr, corresponding 
to the temperature at which the Universe becomes radiation dominated. 

We will consider throughout the case where the reheating temperature from inflaton 
decay satisfies the thermal gravitino constraint [21], Tr ^ 10 8 GeV for 7773/2 = 100 GeV 
( Tr < 10 9 GeV for 777,3/2 = 1 TeV) [22]; in the following we will consider the constraint 
for 7773/2 = 100 GeV. This constraint is not necessary if the curvaton dominates the 
energy density for long enough that it dilutes the thermal gravitinos due to inflaton 
reheating. However, we will be considering the limiting case where the curvaton just 
dominates the energy density before it decays, so we will impose the thermal gravitino 
constraint for consistency. In this case the expansion rate at the end of IMD will 
satisfy H(T R ) < 10~ 2 GeV (1 GeV) for T R < 10 8 GeV (10 9 GeV), which is small 
compared with the value H osc « rriN ^ rn s ~ 100 GeV at which RH sneutrino coherent 
oscillations begin. Here is the RH sneutrino mass, m s is the gravity- mediated soft 
SUSY breaking scalar mass and H(Tr) = k^T^/M is the expansion rate at the onset of 
radiation domination, where k? = (Tc 2 g(T) /9Q) 1 / 2 and g(T) is the number of degrees of 
freedom in thermal equilibrum at T [26 . (In the following we will consider g(T) ~ 200, 
corresponding to all MSSM fields in equilibrium, such that kr ~ 4.7.) Thus we will 
consider the RH sneutrino oscillations to begin during IMD. 

The superpotential of the mass eigenstate RH neutrinos, Ni, is given by 

Mm N 2 

W v = \ v ijNiH u Lj H ^— ^- , (1) 

where H u and Li are the MSSM Higgs and lepton doublets respectively and i is a 
generation index. This leads to Majorana neutrino masses via the see-saw mechanism 
|14j . In addition, we generally expect Planck-scale suppressed non-renormalizable 
superpotential terms involving the RH neutrino superfields of the form 

n!M n-3 ^ > 

In [7] it was shown that such terms must be eliminated for all 77 < 11 if the initial am- 
plitude of the Majorana RH sneutrino curvaton is to be large enough for the coherently 
oscillating curvaton to dominate the energy density before it decays. Otherwise the 



contribution of Planck-suppressed terms to the effective mass squared of the RH sneu- 
trino will be greater than H 2 , thus driving the RH sneutrino amplitude to small values. 
This elimination of the Planck-suppressed NRT requires an explanation in terms of a 
symmetry or other mechanism. One possibility is an R-symmetry. The superpotential 
Eq. is R-symmetric under a simple extension of the R-symmetry of the MSSM 
under which R(Q, L,u c , d c , e c ) = 1 and R(H U , H d ) = [T7|. By assigning R(N) = 1, 
all the terms in Eq. ([TJ are permitted by the R-symmetry whilst all superpotential 
terms in Eq. (0) with n ^ 2 are excluded. The fiH u Hd term, essential for satisfac- 
tory electroweak symmetry breaking ^7], breaks this R-symmetry to R-parity. The 
Giudice-Masiero mechanism [IS] initially eliminates the /i-term via an R-symmetry. 
The //-term is then generated by R-symmetry breaking due to a non-minimal kinetic 
term in the Kahler potential of the form z^H u Hd, coupling the SUSY breaking hidden 
sector superfield z to the Higgs doublets. Since R(z) = 0, this term is permitted by the 
R-symmetry. The F-term of the z superfield breaks both SUSY and the R-symmetry, 
resulting in a //-term of the order of the electroweak scale in the superpotential of the 
low energy effective theory [TSj . The important feature of this mechanism is that it 
only generates R-symmetry breaking terms in the low-energy effective superpotential 
which have total R-charge equal to zero. As a result no superpotential terms of the 
form N n (n ^ 2) will arise. However, this mechanism does require that some interac- 
tions between the hidden sector and observable sector superfields are eliminated by a 
mechanism other than R-symmetry, since R-symmetry permits a superpotential term 
of the form zN 2 , which would typically lead to a very large RH neutrino mass due to 
the z expectation value ~ M. 

The scalar potential of the mass eigenstate RH sneutrinos after inflation is given 

W = (^)W^(^ + 4 (3) 

where we have suppressed the generation indices. Here m 2 N = m 2 + Mjy + m,gjj, where 
m eff = < 2 > is the effective mass squared term due to the interaction of the 
RH sneutrino with thermal background MSSM fields, <fi j2ZJ- The latter will be shown 
to play no role in the evolution of the RH sneutrino curvaton. Since there is little 
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motivation to study RH neutrino masses less than the order of 100 GeV, for simplicity 
we will restrict attention to the case Mjy ^ m s in the following and therefore consider 
rriN ~ Mat throughout. In fact, we will show that the allowed values of Mn are 
generally larger than 0(100) GeV in models where the curvaton condensate decays 
before the electroweak phase transition, which are the most promising with respect to 
baryogenesis. 

The mass squared terms and possibly also B-terms gain corrections during IMD 
due to non-zero inflaton F-terms I2H1 ED]- (During inflation such terms must be 
effectively zero in order to have a scale-invariant density perturbation spectrum, which 
is natural in D-term inflation [TT].) We will consider a negative order H 2 correction 
after inflation with c ~ 1. This is because the evolution of the RH sneutrino under 
the influence of a positive mass squared term in the case c ~ — 1 will suppress the 
initial amplitude of the RH sneutrino oscillation relative to its amplitude at the end 
of inflation. As we show later, this would require a large value of N during inflation 
in order for the curvaton to dominate the energy density before it decays, which 
would in turn require a large value of Hi in order to have sufficiently large curvaton 
energy density perturbations. However, this large value of Hi will generally result in 
excessively large inflaton energy density perturbations and gravitational waves. 

The RH sneutrino oscillations begin at H osc m Mn/c 1 ^ 2 « Mat. The energy density 
of the RH sneutrino condensate at a > or is then [7j 



where an is the scale factor at the onset of radiation domination, p r is the radiation 
energy density, p w 3H 2 M 2 is the energy density due to the coherently oscillating 
inflaton field during IMD and = Mf^N 2 sc /2 is the initial energy density in the RH 
sneutrino oscillations. RH sneutrino domination therefore occurs at jZj 



where H oc a -3 / 2 during IMD. In this expression we have defined curvaton domination 
to correspond to p^ > 4p r , in order to ensure complete domination of the energy 
density by the curvaton. 






(5) 
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We consider first the case of a single neutrino generation, which will later cor- 
respond to the RH sneutrino mass eigenstate with the weakest Yukawa couplings in 
Eq. ((TJ. The RH sneutrino decay rate is given by 



XlM N = m u M 2 N 
Ait 4ttv 2 



where the neutrino mass is given by the see-saw expression 



(6) 



and < H u >= v/2. The decay temperature of the RH sneutrino condensate is then 

^ ( \IM N M \ 1/2 _ ( m u M 2 N M \ 1/2 

Td ~{~^h~) = { Ank T v 2 ) ■ (8) 

Therefore 

Mn „ 160 firaV' 4 f_^_) (lo^evy 2 1 — « — \ GeV (9) 

V 200 / V100GeV/V m v / V 100 GeV J KJ 

From Eq. Q we see that if > 100 GeV and m v is not much smaller than 10~ 4 eV 
then Td > 100 GeV. In this case there is a possibility of generating the baryon asym- 
metry either via electroweak baryogenesis [2H] or via leptogenesis by RH sneutrino 
decay [T3| . 

The condition for the RH sneutrino condensate to dominate the energy density 
before it decays, Hd om > Td, is given by 

Therefore 

N ~ > 4.0 X I0-/J ; = f^V' 4 ( l -^)" 2 (^f . (11) 



M ' ' ' V10- 4 eVy \ T R J 

Thus if the mass of the lightest neutrino is not very much smaller than 10 -4 eV then 
N osc cannot be very much smaller than M. 

A value of N larger than M implies that pn ~ —H 2 N 2 pa —p, leading to a break- 
down of the assumption that the energy density is dominated by inflaton oscillations 



after inflation. In addition, chaotic initial conditions with all initial gradient energy 
densities < M 4 would imply that no field can have a value at the onset of inflation sig- 
nificantly larger than M [9 . We will therefore require that N osc < M. From Eq. (|11|) 
a range of N osc compatible with curvaton domination then exists only if 

A second condition that the RH sneutrino curvaton must satisfy is that the curva- 
ton energy density perturbations are of the correct magnitude to account for the CMB 
temperature fluctuations. The curvaton energy density perturbation is given by [3] 

where Nj is the amplitude of N during inflation and SN w Hi/27r is the quantum 
fluctuation of the RH sneutrino. The observed perturbation is given by 1 5 P ~ 1.4 x 
10~ 4 3J. If N does not evolve from the end of inflation then Nj « N osc > 10~ 3 M, 
which would imply that Hi w 7rS p N osc > 10 12 GeV. However, it was observed in [7] 
that the RH sneutrino field will roll to larger values if the expected negative order H 2 
mass squared correction arises after the end of inflation. In this case N osc 3> Ni is 
possible, allowing smaller values of Hi to be compatible with a RH sneutrino curvaton. 
During IMD, H oc a~ 3 ^ 2 . In this case the evolution of the RH sneutrino due to a —cH 2 
mass squared term is given by iV oc a 7 , where |7j 



1 

Thus 



3 9 , 
-2 + V4 +4c 



(14) 



Nosc = (OoscV = (HW 1 '* ( 15 ) 

Ni \aj \H os J ' 1 ' 

where a osc , a e are the scale factors at the onset of curvaton oscillations and the end of 

inflation respectively. The amplitudes of iV and SN due to a mass squared correction 

evolve in the same way, therefore 5N/N remains constant and Eq. ()13|) is unaltered 

1 More precisely, the power spectrum of the curvaton fluctuation is Vsn/n = Hi/2ttNi = 3V^ 2 /2, 
where the observed value of curvature spectrum is Vl' 2 = 4.8 x 1CT 5 0. 
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[Zj. The curvaton domination condition Eq. (fTTj) combined with Eq. (pH?|) and Eq. (|TH|) 
then implies that the lower bound on Hi is given by 

#1+27/3 > 4X) x 10 _s faSpH^M . (16) 

The lower bound depends on c. If c = 1, for example, then 7 = 0.5 and with 5 P = 
1.4 x 10~ 4 the lower bound on the expansion rate during inflation becomes 

jj > q o ln9 ( Y /W f l ° 8 GeV \ V8 ( M N \V 8 ZlOOGeVX 3 / 8 

(17) 

Smaller values of Hi are possible with larger values of c. For example, with c = 2.5, 
corresponding to 7 = 1, the lower bound becomes 

m„ \ 3 / 20 /i0 8 GeV\ 3/1 °/ Mn \ 7 ' w ( 100 GeV\ 3/10 



^- 2 - 4xl ° 8 llF^vJ l^^J llOO^evJ GeV - 

(18) 

For the case of a positive H 2 correction to the RH sneutrino mass squared with 
c w -1, we see from Eq. (JTIJ) that Re(y) = -3/4 and so N osc /Ni = (H os J Hi) 1 / 2 » 
(M N /Hi) 1/2 . With JVj related to if/ by Eq. ffT3j). iiV is therefore given by 

7T 2 (5 2 iV 2 

*' = -ikf ■ (19) 

The lower bound on N osc , Eq. (|11|) . then implies that 

V10- 4 eV7 V T « / 

if/ is therefore much larger than the upper bound from gravitational waves produced 
during inflation, Hi < 10~ 5 M |27)j . 

One important application of a SUSY curvaton is to D-term hybrid inflation jllj . 
It has been shown that the contribution of U{1)fi strings to the CMB, assuming 
naturally large values of the order of 1 for the superpotential coupling, is inconsistent 
with observations in the case where the energy density perturbations are explained by 
conventional inflaton perturbations, in order to evade this problem, the expansion 
rate during inflation must satisfy Hi < 4 x 10 12 g GeV |12j . where g is the Fayet- 
Iliopoulos gauge coupling [TT] . In this case a new source of energy density perturbations 
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is required, since conventional inflaton energy density perturbations are too small. 
From Eq. (|17|). requiring that Hi < 4x 10 12 g GeV implies that a RH sneutrino curvaton 
with c = 1 can evade the D-term inflation cosmic string problem and provide the energy 
density perturbations if 

whilst for c = 2.5, 

M^,0 Xl 0V-(^) 3/ "(^) S/ '(^)- G eV. <*> 

3 Thermal Constraints on Majorana Right-Handed 
Sneutrino Curvatons 

In this section we discuss the constraints coming from the requirement that the Ma- 
jorana RH sneutrino condensate is not thermalized before the condensate decays. We 
also comment on the possibility that an effective RH sneutrino mass due to the inter- 
action of the curvaton RH sneutrino with the thermal background could play a role in 
the evolution of the condensate j27j . 

In [Z| it was observed that a Majorana RH sneutrino condensate will be thermalized 
before it decays if the Yukawa couplings of the mass eigenstate RH neutrinos are all 
of the same order of magnitude. This would be expected if there was no correlation 
between the RH neutrino mass matrix and the Yukawa coupling matrix X U ij, in the 
sense that the basis in which the RH neutrino mass matrix is diagonalized is unrelated 
to the basis in which the Yukawa coupling matrix is diagonalized, since in this case all 
the mass eigenstate RH neutrino superfields would have similar Yukawa couplings to 
the lepton doublets. The thermalization rate would then be determined by the largest 
Yukawa coupling, corresponding to the neutrino mass indicated by atmospheric solar 
neutrino oscillations, m v « 0.1 eV. 

However, it is possible for one of the RH neutrino mass eigenstates to have much 
weaker Yukawa couplings to the MSSM fields than those which are responsible for 
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the mass splittings and mixing angles observed in solar and atmospheric neutrino 
oscillations JE] • in this case there will be a hierarchical light neutrino mass spectrum 
m Vl <C m V2 <C Tn U3 , where m Vi (i = 1,2,3) are the light Majorana neutrino mass 
eigenvalues ^5] • It is then possible for coherent oscillations of the weakly coupled RH 
sneutrino to evade thermalization. Since coherent oscillations of the more strongly 
coupled RH sneutrinos will be thermalized by scattering from the thermal background 
MSSM fields, the weakly coupled RH sneutrino will naturally play the role of the 
curvaton. 

For simplicity, we will consider the Yukawa couplings of the weakly coupled RH 
sneutrino mass eigenstate N% to lepton doublets Li to be all of the same order of 
magnitude, A„ u ~ X u . X u will then determine the lightest neutrino mass m ui via the 
see-saw mechanism with N\. It will also determine the thermalization rate of the RH 
sneutrino condensate. We will drop the generation subscript and write N\ as N in the 
following. 

The condition for the RH sneutrino condensate to evade thermalization is that the 
scattering rate of the most strongly coupled thermal background MSSM fields from 
the zero-momentum condensate RH sneutrinos is less than the expansion rate H. The 
strongest scatterings will be from top quarks/squarks in the thermal background via 
Higgs exchange. The scattering cross-section is given by [7] 

a u a t oi v a t 



E 2 M M N T 



(23) 



where a Ut % = X 2 , t /47r, At is the top quark Yukawa coupling and E^ M = txi^Et is the 
centre of mass energy of the scattering particles, where Et ~ T is the energy of the 
thermal background particle. The scattering rate from relativistic top quarks/squarks 
is then T sc = na sc v = na sc where n ~ g t T 3 /7r 2 is the number density of top quarks 
and squarks (g t ~ 10 is the number of degrees of freedom in top quarks and squarks 
including spin and colour) and v = 1 for relativistic thermal particles. Thus 

r « = ^M^ < 24) 

The important feature of this expression is that the scattering rate is proportional to 
T 2 , as is H during radiation domination. Therefore the condition for the condensate to 
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evade thermalization by scattering during radiation domination, T sc < H = kxT 2 /M, 
is independent of the temperature, 

a; : . (25) 

Using the relation between the Yukawa coupling X u and the lightest neutrino mass 



m 



V\ ■ 



Eq. (J7J), then implies that 

. 4vr 3 fc T v 2 o /10\ /0.1\ / v \ 2 

m Ul < -— » 2 x 10~ 3 — — 7T — T eV , 26 

1 g t a t M \g t J \a t J U00 GeV J v ; 

where we have used &t = 4.7. From this we conclude that for m Vl ^ eV it should 
be possible for the Majorana RH sneutrino curvaton to evade thermalization during 
the radiation dominated era, T < Tr. 

This condition is also sufficient to ensure that the curvaton is not thermalized by 
the background radiation coming from inflaton decays during the IMD era. To see this, 
note that the expansion rate is related to the background radiation temperature during 
IMD by T = k r (M Pl HTlf ,A where k r = (9/5tt 3 ^(T)) 1 / 8 [2E|. Thus H oc T 4 whilst 
T sc is proportional to T 2 as before. Therefore as H and T increase, the scattering 
rate becomes small relative to the expansion rate. Thus the largest scattering rate 
relative to the expansion rate will occur at the lowest temperature during IMD, Tr. 
However, we know that T sc < H at Tr if Eq. ()26|) is satisfied. Therefore if Eq. (|26j) 
is satisfied, the RH sneutrino condensate will evade thermalization during both the 
radiation dominated and IMD era. 

A second thermal effect is the effective mass term generated by the interaction of 
the RH sneutrinos with the thermal background MSSM fields [21] ■ This gives rise 
to a mass squared term of the form X 2 < <p 2 >m \ 2 T 2 , where is a H u or Lj field 
in the thermal background and < ... > denotes the thermal average. Requiring that 
this is small compared with implies that \ v T max < Mjv, where T max is the largest 
background temperature experienced by the coherently oscillating RH sneutrino. This 
corresponds to the temperature at the onset of RH sneutrino oscillations, T max = 
k r (M Pl H osc T 2 ) 1/4 w k r (M Pl m N T 2 ) 1/4 , where k r w 0.4 for g(T) = 200. Therefore 
il4 /4 7 f M N \ 3 / 4 /i0 8 GeV\ 1/2 
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This translates into an upper bound on the lightest neutrino mass 

m - 2Xl ° (iMGiv) (lOOGiv) HH^' (28) 

Thus m Vl < 10~ 3 eV both allows the Majorana RH sneutrino condensate to evade 
thermalization and ensures that thermal effective masses will play no role in the evo- 
lution of the RH sneutrino curvaton. The existence of a very light neutrino eigenstate 
may be regarded as a prediction of the Majorana RH sneutrino curvaton model. 

4 Lepton Asymmetry Constraints 

If the RH sneutrino condensate decays at Td > T ew then we must ensure that any lepton 
asymmetry produced by the decay of the condensate, which would be converted into 
a baryon asymmetry of equal order of magnitude by sphaleron processes (ub/s = 
—ni initial /2s [2E]), is not larger than the observed baryon asymmetry of the Universe. 

In general, a coherently oscillating complex scalar field in SUSY will acquire a global 
charge asymmetry via the Affleck-Dine (AD) mechanism [221 • Soft SUSY breaking A- 
and B-terms in the scalar potential [Uj will cause the real and imaginary parts of the 
complex scalar field to oscillate out of phase, such that the scalar field descibes an el- 
lipse in the complex plane corresponding to a globally charged scalar field condensate 
|221 ITU] . However, it is not possible to generate the baryon asymmetry via AD baryo- 
genesis or leptogenesis from the decay of a curvaton, since in this case there would 
be large baryon isocurvature perturbations due to the fact that in the AD mechanism 
both the real and imaginary parts of the curvaton field act to determine the baryon 
or lepton asymmetry [3TJ. 

The global U(l) symmetry in the case of a RH sneutrino condensate is lepton 
number L whilst the B-terms driving the Affleck-Dine mechanism are those associated 
with the RH neutrino mass terms in Eq. (|Tjl. The B-term in Eq. (jSJ) splits the masses 
of the real and imaginary parts of the RH sneutrino field. For H 2 <C the scalar 
potential is given by 

N 2 N 2 

V{N) = {M 2 N + BM N ) + [m 2 n — BM N ) -f , (29) 
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where N = (N\ + iN 2 )/\/2 and we have taken the RH sneutrino mass m^ in Eq. (J3j) 
to be equal to Mat- (We may choose the phase of N such that B is real and positive.) 
The solutions for coherent oscillations of the real and imaginary parts of the RH 
sneutrino field are then N\ = A\ cos (mit) and N 2 = A 2 cos [m 2 t + 5), where mi = 
(Ad% + BM N ) 1/2 , m 2 = (M 2 N - BM N ) l/2 and A h 2 oc a~ 3 / 2 . We have included a phase 
shift S, since this would be generally expected to arise during the initial transition from 
Hubble damped evolution to coherent oscillations. We will see that S has a negligible 
effect on the net lepton asymmetry. Since the B-term in the scalar field equation for 
N becomes dynamically significant, BMn Z H 2 , only once H « Mat, corresponding 
to the onset of coherent oscillations, there is no reason to expect N to be along the 
real direction at the onset of RH sneutrino oscillations. Thus we expect the phase of 
N initially to be random and of the order of 1, such that A\ and A 2 will be of the 
same order of magnitude. 

The lepton asymmetry density in the RH sneutrino condensate, til cond, is given by 

n L ^ = % (N^N - N^N) = (N 2 N, - jViJV 2 ) . (30) 

Therefore 

n L cond = A\A 2 (mi sin(mit) cos(m 2 t + 5) — m 2 sm(m 2 t + 5) cos(mit)) . (31) 

For example, for the case 5 = 0, if we consider 3> B then mi ~ Mjy + 5m/2, 
m 2 ~ Mjv — dm/2, where 5m w B. Then 

(OTYl \ 
sin (5mt) + — sin (2mt) I . (32) 

Thus, for m, 8m 3> H(T<i), there will be rapid oscillations of the lepton asymmetry 
between large positive and negative lepton number over the time scale of condensate 
decay, I^ 1 w H^{T d ). 

We next calculate the net lepton asymmetry produced by the decay of the RH 
sneutrino condensate. Although the lepton asymmetry in the condensate is time 
dependent, once lepton number is transferred to the thermal background it will be 
conserved. (We neglect the effect sphaleron processing for now, which will not affect 
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the magnitude of the lepton asymmetry.) This follows since lepton number violating 
scattering processes due to RH neutrino exchange are out of equilibrium at tempera- 
tures below Tf 1.2 x 10 11 (1 eV/m u ) 2 GeV [32|- Therefore the rate of increase of the 
thermal background lepton number density, n^, from RH sneutrino condensate decay 
is given by 

1 rl ( n^n t\ 

= T d e~ d n L cond , (33) 



a 3 dt 

where we have included a factor e~ Tdt in order to account for the decay of the RH 
sneutrino condensate and the scale factor terms account for the expansion of the 
Universe. Thus 

riL = T d AiA 2 / e~ Vdt [(mi sin(mii) cos(m 2 t) — m 2 sin(m 2 £) cos(rait)) cos S 
Jo 

— (mi sin(mit) sin(m2t) + m 2 cos(mit) cos(m2t)) sin 5} dt . (34) 

The integrals may be calculated exactly, with the result 

T d A 1 A 2 (ml - ml) [T 2 d + m\ + m£] cos 5 



T 2 , + (mi - m 2 ) + (mi + m 2 ) 
r^iv4 2 (m 2 r^ + 3mjm 2 + m|) sin 5 



T 2 , + (mi - m 2 ) + (mi + m 2 



(35) 



For small T 2 d compared with (mi — m 2 ) 2 , and assuming that the natural value of 5 is 
not very close to ir/2, the term proportional to sin 5 is negligible and Eq. (j33j) becomes 



r d AiA 2 (m 2 + m 2 ) cos 5 _ TjA^Mn cos 6 
HL ^ (ml - ml) " B ' 



We wish to express this in the form n^js. The energy density in the RH sneutrino 
condensate is given by 

•m\A\ m\A\ 

P = pi+p2 = + ■ (37) 

The entropy density at T d is related to p by s = 4p/3T d j2Sj, where we assume the 
radiation energy density from RH sneutrino decay at T d equals p. Thus the lepton 
number to entropy density is given by 

«l _ 3T d T d M N cos5 A X A 2 

s ~ 2B (m 2 A 2 + m 2 A 2 ) ' 1 j 
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With A\ ~ A 2 , ml pa ml pa M%, cos 5 pa 1 and < M^r, this becomes 

n L T d T d 



(39) 



Using T d = kxT^/M, the lepton asymmetry from B-term AD leptogenesis is therefore 
given by 



s BM N M 
Requiring that ul/s < 1CT 10 imposes the constraint 



(40) 



With T d given by Eq. (??), this implies an upper bound on M 



N 



This may also be expressed as an upper bound on T d , 

< 8 . 6 x 10 . (i^f ( _4_ ) 1/2 GeV . (43) 
V m ui J VIOOGeV/ v ; 

Thus if m ui is not very much smaller than the mass of the other neutrino eigenstates 
then the lepton asymmetry constraint Eq. (|42|) favours a relatively light RH neutrino 
mass, rriN ^ 10 5 GeV, whilst Eq. (J43|) implies that the decay temperature is relatively 
low, T d < 10 5 GeV. 

The analysis of the B-term AD mechanism given above is similar to that recently 
given in j^j in the more general context of AD leptogenesis from a RH sneutrino 
condensate. Some of the analytical results we have obtained here should be useful 
in this more general context, in particular the expression for the lepton asymmetry 
resulting from decay of the condensate, Eq. (|3*3]l. From Eq. (J33j) we see that the 
maximum asymmetry will be obtained for T d pa \ mi — 777.2 1 ~ B, in agreement with [21J. 
However, it is worth emphasizing that a significant asymmetry can be still obtained 
even if T d -C B, as is clear from Eq. (|39|) . 

There are two obvious possibilities for the origin of the baryon asymmetry in the 
Majorana RH sneutrino curvaton model: electroweak baryogenesis [231 or RH sneu- 
trino leptogenesis based on conventional leptogenesis via CP violating decays of the 
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condensate RH sneutrinos [13J. Both of these require that the curvaton condensate 
decays before the electroweak phase transition. In the case where the energy density 
is dominated by the RH sneutrino condensate, it was shown in that the lepton 
asymmetry from CP-violating decays of the condensate RH sneutrinos is given by 

- * 0.7 x 10" 10 ( J d ) ( m " ) 6 eff (44) 
s V10 6 GeW V0.05 eV J n V ; 

where 8 e ff < 1 is the CP- violating phase. Therefore RH sneutrino leptogenesis requires 
that Td m l&GeV if 5 e ff ~ 1. From the lepton asymmetry constraint, Eq. (|41j). this 
would require that M N > 10 8 GeV. Eq. in turn requires that m Ul < 10~ 8 eV. 
Therefore RH sneutrino leptogenesis is only possible if the lightest neutrino mass is 
extremely small compared with the mass scale suggested by solar and atmospheric 
neutrino oscillations. 



5 Consequences of the Constraints 

We next consider the consequences of the simultaneous application of the constraints 
derived above for the parameter space of the Majorana RH sneutrino curvaton model. 
We will consider the allowed region of the (Mat, Tr) parameter space for different values 
of m Vl satisfying the thermalization upper bound m Ul < 10~ 3 eV and for various lower 
bounds on T^. 

The constraints we apply are: (i) the upper bound on from curvaton dominance 
combined with N osc ^ M, Eq. (|12|) . (ii) the lower bound on Mjv from requiring that the 
RH sneutrino decay temperature in Eq. Q is greater than a given value and (iii) the 
upper bound on from the AD leptogenesis constraint, Eq. (|4"2jl . In addition, we will 
consider the upper bound on Mjv for which a RH sneutrino curvaton is compatible with 
a solution of the D-term inflation cosmic string problem for the case c = 1, Eq. (|21|1. 
(We consider g — 1 for the U(1)fi gauge coupling throughout.) 

In Figure 1 we show the constraints for the case m Vl = 10~ 3 eV (the largest value 
consistent with thermalization) and lower bound on equal to 100 GeV, correspond- 
ing to the value for which electroweak baryogenesis is possible. The case m Ul & 
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Figure 1: Right-handed neutrino mass vs. reheating temperature for m vi = 1CP 3 eV and T d > 
10 2 GeV. The solid line is the lower bound on M N from T d > 10 2 GeV. The long-dashed line is the 
upper bound from B-term AD leptogenesis. The dot-dashed line is the upper bound from curvaton 
dominance. These define the allowed region, shown hatched. The short dashed line is the D-term 
inflation upper bound on Mm for c = 1 . 
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Figure 2: Right-handed neutrino mass vs. reheating temperature for m Vl =10 4 cV and > 
10 3 GeV. 
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,'U.re 3: Right-handed neutrino mass vs. reheating temperature for m Vl — 10 9 eV and Td 
GeV. The D-term inflation upper bound for the case c = 2.5 is also shown (dotted line). 
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1CT 3 eV might be expected if m ui is not very much smaller than m V2 3 ~ 0.01 — 0.1 eV. 
From Figure 1 the allowed range of Mm is constrained to be between 47 GeV and 
2.3 x 10 4 GeV. The allowed range of Mn becomes narrower for Tr < 10 6 GeV, where 
the curvaton dominance constraint becomes more important than the AD leptogen- 
esis constraint. Tr cannot be smaller than 2.7 TeV. It is interesting and perhaps 
significant that the allowed range of M^r contains the 100 GeV — 1 TeV mass range, 
a natural mass range of particle physics. 

In Figure 2 we consider a more restrictive possibility with m Vl = 10~ 4 eV and 
Td > 1 TeV. The allowed range of Mjv is significantly narrower in this case compared 
with Figure 1, with M N between 1.4 x 10 3 GeV and 1.3 x 10 5 GeV. In general, 
increasing the lower bound on Td increases the lower bound on Mjy (oc Td, Eq. Q) 
whilst having no effect on the upper bound. Decreasing m vi raises both the upper 
bound (oc m~ 3//4 , Eq. and, to a lesser extent, the lower bound (oc m" 1 / 2 Eq. 0). 

In Figure 3 we show the case where the lower bound on Td is 10 6 GeV, compatible 
with RH sneutrino leptogenesis. In this case the lower bound on Mjv is 5.0 x 10 s GeV. 
We have shown the allowed parameter space for the case where m Vl = 10~ 9 eV, which 
is just compatible with RH sneutrino leptogenesis, as illustrated by the very small 
region of the parameter space compatible with the constraints. For m Ul > 10~ 8 eV 
RH sneutrino leptogenesis is ruled out by the AD leptogenesis constraint. The case 
Td = 10 6 GeV is not compatible with a curvaton solution of the D-term inflation 
cosmic string problem when c = 1. However, larger values of c (but still of the order 
of 1) will allow a solution, as shown by the upper bound for the case c = 2.5. 

We conclude that unless m Ul is very small compared with the scale of neutrino 
masses responsible for solar and atmospheric neutrino oscillations, the Majorana RH 
sneutrino curvaton favours relatively low values of and T^. The region of the 
parameter space compatible with a Majorana RH sneutrino curvaton is generally 
consistent with a curvaton solution of the D-term inflation cosmic string problem. 
RH sneutrino leptogenesis is ruled out as a source of the baryon asymmetry unless 
m u 1 ^ 10~ 8 eV, thus favouring electroweak baryogenesis as the origin of the baryon 
asymmetry. However, it may be significant that RH sneutrino leptogenesis is not 
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entirely ruled out as a source of the baryon asymmetry. 

6 Conclusions 

We have considered the conditions under which a Majorana RH sneutrino could serve 
as a curvaton. One condition is the absence of Planck-scale suppressed corrections to 
the RH sneutrino scalar potential to a high-order. As the nature of such corrections 
requires speculation about the nature of the underlying full theory, we have considered 
as one possibility the case where such corrections are absent. This might be explained 
by an R-symmetry, although this requires an additional suppression of interaction 
terms between the observable sector and SUSY breaking hidden sector superfields 
in the superpotential. A second condition is that the RH sneutrino condensate is 
not thermalized by scattering from the thermal background particles. This requires 
a model of neutrino masses in which one of the RH neutrino mass eigenstates has 
significantly weaker Yukawa couplings to the MSSM fields than the other RH neutrinos. 
The lightest neutrino mass eigenstate must then satisfy m ui < 1CT 3 eV, which may 
be regarded as a prediction of the scenario. The Majorana RH sneutrino curvaton 
scenario can therefore be ruled out experimentally by the observation of a degenerate 
neutrino mass spectrum, m vi m V2 m„ 3 . 

We have shown that a lepton asymmetry will be induced in the RH sneutrino con- 
densate via the AD mechanism associated with SUSY breaking B-terms. Requiring 
that the baryon asymmetry coming from sphaleron conversion of this lepton asym- 
metry is smaller than the observed baryon asymmetry imposes an upper bound on 
M N . Unless m vi is very small compared with the mass scale of solar and atmospheric 
neutrino oscillations, relatively small values of M N ( < 10 5 GeV) are required for con- 
sistency with the AD leptogenesis constraint. Some of the analytical results we have 
obtained here should be useful in the more general analysis of AD leptogenesis from a 
RH sneutrino condensate. 

The AD leptogenesis constraint, combined with a lower bound on the RH sneutrino 
condensate decay temperature and the requirement that the RH sneutrino curvaton 



24 



condensate dominates the energy density before it decays, restricts the region of the 
(Mn,Tr) parameter space consistent with a RH sneutrino curvaton. The range of 
Mn consistent with a RH sneutrino curvaton is typically smaller than that considered 
in the conventional see-saw mechanism where A„ is assumed to be of the order of 
the charged lepton Yukawa couplings, which requires « 10 9 GeV. However, the 
allowed range of for m Vl ~ eV and Td > 100 GeV is consistent with the 
100 GeV — 1 TeV mass scale, a natural mass scale in particle physics models. 

An important motivation for a SUSY curvaton is D-term hybrid inflation. We 
have shown that the region of the parameter space allowed by the constraints on 
the Majorana RH sneutrino curvaton is generally consistent with a solution of the 
D-term inflation cosmic string problem. D-term inflation naturally accounts for the 
absence of order H 2 corrections during inflation in both the inflaton and curvaton 
scalar potentials. Therefore a combined D-term inflation/Majorana RH sneutrino 
curvaton model could provide a completely consistent model of SUSY inflation. 

The only obvious candidates for the origin of the baryon asymmetry in the Majo- 
rana RH sneutrino curvaton scenario are electroweak baryogenesis and RH sneutrino 
leptogenesis, AD leptogenesis being generally ruled out for a curvaton |3~T] . We find 
that RH sneutrino leptogenesis via CP violating decays of the RH sneutrinos is con- 
sistent with the Majorana RH sneutrino curvaton only if the lightest neutrino mass 
eigenstate satisfies m Vx < 10~ 8 eV. Therefore RH sneutrino leptogenesis appears to 
be disfavoured, although it may be significant that it is not completely ruled out. 
This suggests that electroweak baryogenesis is the most likely source of the baryon 
asymmetry in this scenario. 
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